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We study the effect of electron-electron interactions in the electronic properties of a biased 
graphene bilayer. This system is a semiconductor with conduction and valence bands character- 
ized by an unusual "mexican-hat" dispersion. We focus on the metallic regime where the chemical 
potential lies in the "mexican-hat" in the conduction band, leading to a topologically non-trivial 
Fermi surface in the shape of a ring. We show that due to the unusual topology of the Fermi surface 
electron-electron interactions are greatly enhanced. We show that the ferromagnetic instability can 
occur provided a low density of carriers. We compute the electronic polarization function in the 
random phase approximation and show that, while at low energies the system behaves as a Fermi 
liquid (albeit with peculiar Friedel oscillations), at high frequencies it shows a highly anomalous 
response when compare to ordinary metals. 

PACS numbers: 73.63.-b, 71.70.Di, 73.43.Cd 
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I. INTRODUCTION 

The discovery of graphene^, and the reahzation that 
it presents very unusual electronic properties^, has gen- 
erated an enormous amount of interest in the condensed 
matter community. The possibihty of creating electronic 
devices consisting of only a few atomic layers can open 
doors for a carbon-based micro-electronics. In this con- 
text, bilayer systems'^ play a distinguished role for being 
the minimal system to produce a gap in the electronic 
spectrum due to an electric field effect^'^. Recent ex- 
periments show that the electronic gap and the chemi- 
cal potential can be tuned independently of each other 
and the band structure can be well described by a tight- 
binding model corrected by charging effects^. The elec- 
tronic gap in these systems has been observed recently 
in angle resolved photoemission (ARPES)^ in epitaxially 
grown graphene films on SiC crystal surfaces^. 

Electron-electron interactions are usually neglected in 
single layer graphene since they do not seem to play 
an important role in the transport measurements in low 
magnetic fields, e. g., in the measurements of the anoma- 
lous quantum Hall effect^. At high magnetic fields, how- 
ever, the electronic kinetic energy becomes quenched by 
the creation of Landau levels^*^ and Coulomb interac- 
tions become important"'^^'^^. An alternative theory pre- 
dicts that Coulomb interaction becomes relevant even 
at arbitrary small magnetic fields"'^^. Also, the minimal 
conductivity^^ and collective excitations"'^'* depend on the 
electron-electron interaction. 

In single graphene sheets, the unscreened electron- 
electron interaction is - due to the vanishing density of 
state at the Dirac point - marginally irrelevant in the 
renormalization group sense. This leads to a behav- 
ior with the electronic self-energy of the form, ImE cx 
uj,^^ which is reminiscent of so-called marginal Fermi 



liquid behavior^^. The fact that even the unscreened 
Coulomb interactions are marginally irrelevant in sin- 
gle layer graphene suppresses the presence of magnetic 
phases in this system. In fact, ferromagnetism can only 
be found in strong coupling, when graphene's fine struc- 
ture constant, ac = e^/ie-ohvp) (e is electric charge, eg 
is the dielectric constant, and vp ^s, the Fermi-Dirac ve- 
locity), is larger than a critical value ac- It was found 
that ac ~ 5.3 for clean graphene {uc ~ 3.8 for disor- 
dered graphene) For tQ — 1 we have ac ~ 2 which is 
smaller than the critical value indicating the absence of 
a ferromagnetic transition in single layer graphene. For 
Coulomb effects in disordered graphene, see also Refs. 
19. 

In this paper, we study the electron-electron interac- 
tion in the biased bilayer system (or unbiased graphene 
bilayer see Ref.^'^). We use the standard tight-binding 
description of the electronic structure^ ^ and treat the 
problem in the metallic regime when the chemical poten- 
tial lies in the conduction band leading to a topologically 
non-trivial ring for the Fermi surface. We also disregard 
inter-band transitions by assuming that the electronic 
gap is sufficiently large so that the valence band is com- 
pletely filled and inert. We show that the ferromagnetic 
instability can occur in this system due to the reduced 
phase space of the ring. We also calculate the polariza- 
tion properties of this system and show that they are 
rather unusual at high frequencies with peculiar Friedel 
oscillations at low frequencies. 

The paper is organized as followed. In section II, we 
introduce the model and discuss the self-energy effects 
due to electron-electron interactions. In section III, we 
present a discussion of the stability of the system towards 
ferromagnetism by introducing a Landau-Ginzburg func- 
tional for the free energy. In section IV, we calculate two- 
particle correlation functions, i.e., the imaginary part of 
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FIG. 1: (Color online) Lattice structure of a bilayer graphene. 
The A and B sublattices are indicated by white and red 
spheres, respectively. 



the polarization and the plasmon dispersion. We close 
with conclusions and future research directions. 



II. THE EFFECTIVE MODEL 

The graphene bilayer consists of two graphene planes 
(labeled 1 and 2), with a honeycomb lattice structure, 
and stacked according to Bernal order (A1-B2), where A 
and B refer to each sublattice of a single graphene (see 
Fig. 1). The system is parameterized by a tight-binding 
model for the 7r-electrons with nearest-neighbor in-plane 
hopping t (« 2.7 eV) and inter-plane hopping t± (~ 0.3 
eV) . The two layers have different electrostatic potentials 
parameterized by V. The electronic band structure of 
the biased bilayer is obtained within the tight-binding 
description^, leading to the dispersion relation (we use 
units such that h = 1 = fc^): 



±± 



iV)^±^el + tl/2 + Vy4±^ti/4+{tl + V^. 
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where et is the dispersion of a single graphene layer. 
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where vp = 3ta/2 where a — 1.4lA is the lattice-spacing. 
The energy dispersion (4) has the shape shown in Fig. 3, 
with its minimum at momentum 



fco = vW(2A), 



(6) 



For a finite density of electrons in the conduction band 
the occupied momentum states are constrained such that 



< A: < fcniax.o- where: 
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"-min.cr — ^ ~ 27m(j , 
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where rio- is the electronic density per unit area for elec- 
trons with spin a — ±1. 

The density of states per unit area p(e) can be written 
as: 



Pie) 



1 



(9) 



with e,jiin = A — q;2/(4A), which has a square root singu- 
larity at the band edge, just like in one-dimensional (ID) 
systems. This ID singularity in the density of states of 
the biased bilayer leads to an strong enhancement in the 
electron scattering either by other electrons or impurities 
(which is beyond the scope of this paper). The Fermi 
energy, ep, can be expressed in terms of the electronic 
density per unit area, as 



(10) 



ek = t[3 + 2cos{akx) -f 4cos {akx/2) cos [afc^ ^3/2]]^/^ (2) The simplest model including the effect of electron- 
electron interactions in this problem is: 



where k — (kx, ky) is the two-dimensional (2D) momen- 
tum. This dispersion relation around the K and K' 
points of the Brillouin zone are shown in Fig. 2. No- 
tice that the electronic spectrum shows a "mexican-hat" 
dispersion with a gap minimum, Ag, given by: 

^,^V[tl/{tl + V^)Y/\ (3) 

that is located at momentum /cq relative to the K point. 
Close to the K point the electronic dispersion can be 
written as: 



H = Y,E{k)cl^Ck., (11) 

k.(T 

where Ck u (c\. ^) annihilates (creates) and electron with 
momentum k and spin a. Notice that only the Coulomb 
interaction between the electrons of the occupied con- 
duction band is considered. The Fourier transform of 
the Coulomb potential is given by: 



E{k) ^A~ak'^ + Xk^ , (4) 
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A. First order corrections to the quasi-particle 
spectrum 

The first order correction to the electronic self-energy 
is given by the diagram in Fig. 4 (a), and it reads: 



Si(fe) 



27re2 



(14) 



FIG. 2: a) The four bands E^'^ (in eV) as function of |k| (in 
A"'^) around the K point with V — t±. b) The "mexican- 
hat" E^'~ (in eV) as function of k (in A^^) around the K 




FIG. 3: "Mexican-hat" dispersion of the biased graphene bi- 
layer close to the K point in the Brillouin zone. The symbols 
are explained in the text. 



The unperturbed electronic Green's function reads: 



(13) 



where w„ = {2n + l)'^/(3 are fermionic Matsubara fre- 
quencies. 




FIG. 4: Evaluated diagrams where the solid line resembles 
the bare electronic propagator and the wiggled line the bare 
Coulomb interaction: a) electronic self-energy, b) particle- 
hole bubble or polarization. 



and considering the zero temperature limit, one gets 



Si(fc) 



27reo 



27reo 



-/l(a, fc/fcmax) , (15) 

1/1 \ P 



+ 32 I rki^. 



where a = fc„iin/fcmax and the Coulomb integral Ii{x) is 
given in appendix A. From this result one obtains the 
renormalized values for A, a and A, defined as Ei{k) — 
E{k) + T,i{k). For k^iin > 0, i.e., the Mexican hat is 
partially filled, the renormalized values read 



Ai — A (^max ^min) j 



ai 



a - 



Ai = A 



e 

4eo 



■\p^ /hi — h'^ 



^max'^inin 



64eo V fcinax^iin 



(16) 

(17) 
(18) 



A consequence of the renormalization of the bare pa- 
rameters is that for densities such that Ai < 0, the 
spectrum becomes unbounded indicating the a possible 
instability of the non-interacting system. For realistic 
values of the parameters (see Section III), the instabil- 
ity occurs at /cmin/fco ~ 3/4, i.e., for densities close to 
Tio- ~ — ^o/C^'"') when the Mexican hat is completely 
filled. We will confine ourselves to the limit of low densi- 
ties so that the expansion of the energy dispersion in Eq. 
(4) up to the quartic term in the electron dispersion is 
valid. For larger densities for which Ai w 0, higher order 
terms in the expansion must be included. For ria > n'^, 
i.e., the Mexican hat is completely filled, the electron- 
electron interaction stabilizes the dispersion and the dip 
in the spectrum reduces. In this case, the parameters of 
the model modify to: 



Ai — A ^max J 

o2 



e 

4eo' 



Cl^l ^ ^max I 



3e2 



(19) 
(20) 
(21) 
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III. EXCHANGE INSTABILITY OF THE 
GROUND STATE 

In this section we examine a possible instability of the 
Fermi ring is unstable towards a ferromagnetic ground 
state. For the instability to occur, the exchange energy 
due to Coulomb interactions has to overcome the loss of 
kinetic energy when the system is polarized. The kinetic 
energy for electrons of spin a has the form: 



(22) 



where A is the area of the bilayer. The exchange energy 
for electrons with spin a can be written as 



^ J2 nik-k'i)i(k|k')i= 



(23) 



k.k' occup. 



Within the approximations described in the introduction, 
i.e., assuming that the wave functions are mostly local- 
ized on one sublattice in one layer, the overlapp factor 
(k|k') « 1. This leads to 



e^A 3 



h ■ 



(24) 



where fcmax.o- and h 



r are defined in Eqs. (7) and (8), 
(see Fig. 3). 

In the paramagnetic phase both spin bands have the 
same number of electrons and hence Ua = n/2. In the 
ferromagnetic phase the system acquires a magnetiza- 
tion density, m = 7i+ — n_, and the electron occupation 
change to: — {n + am)/2. In order to study a possible 
ferromagnetic transition, we derive a Landau-Ginzburg 
functional in powers of m, with the coefficients depend- 
ing on n. The Landau-Ginzburg functional is defined as: 

F[a, A, n, to] = ^[K„{n, to) -|- E^''''{n, m) (25) 

(T 

~ a{n)rn? + b{n)m'^ + c{n)'mP -|- . . . . 

A second order magnetic transition occurs when a(n) be- 
comes negative. If a(n) is always positive, a first order 
transition may occur if b{n) becomes negative. The term 
linear in Ua- in the kinetic energy (22) does not contribute 
to F[a, X, n, to]. 

The coefficient a{n) has two main contributions: one 
comes from the kinetic energy (22), a^{n), and another 
arises from the exchange (24), a°^(n). These contribu- 
tions read: 



pi A 

4eo 



Li +L2 



(26) 
(27) 




FIG. 5: 
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(Color online): Plot of a{n)/A (in units of 
~^]~^), as function of the electronic density n (in 



units of 10 cm~ ). The dashed (green) line gives the bare 
value of a{n)/A, the full (red) line contains self-energy cor- 
rections. 



where the coefficients Li are functions related to the 
Coulomb integral hix) defined in appendix A in the fol- 
lowing way: 



Li 



■ max^ 



fr2 



^ "-max "-n 

97,3 VI 
^""max '^min 



kl(kl - 2mr) 

01,2 p ' 
^'^max "-min 



In the following, we choose the parameters as in section 
II with V — t±. For ac ~ 2 all coefficients are negative 
indicating that the system has a strong tendency towards 
ferromagnetic ordering. In this case the transition is of 
second order. There are, however, values of ac for which 
a(n) is positive, whereas the remaining coefficients are 
negative. In this case the transition is of first order. At 
zero temperature the magnetization density equals the 
density of electrons. 

We have also investigated how the above picture 
changes if one includes the corrections to energy spec- 
trum due to the electron-electron interaction, given by 
Eqs. (17) to (18). The comparison is shown in Fig. 5 
and Fig. 6 where the dashed line shows the coefficients 
without self-energy effects whereas the full line includes 
self-energy effects. The inclusion of self-energy effects 
decreases the tendency toward a ferromagnetic ground- 
state and there is a critical density for which both co- 
efficients, a{n) and 6(n), become positive. This is also 
seen in Fig. 7, where the free energy F{m) is shown for 
various electronic densities. Including self-energy effects 
leads to the saturation of the magnetization at high den- 
sities, i.e., for n = 5 X lO^^cm"^ the magnetization is not 
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IV. POLARIZATION 
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FIG. 6: (Color online): Plot of b{n)/A (in units of 
eV[10^^cm~^]~"^), as function of the electronic density n (in 
units of lO^^cm"^). The dashed (green) line gives the bare 
value, the full (red) line is the result with self-energy correc- 
tions. 
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FIG. 7: (Color online): Plot of F{m)/A (in units of 
eV[10^^cm~^]), as function of the magnetic density m (in 
units of lO^^cm"^) for various electronic densities n includ- 
ing self-energy corrections. The dotted-dashed (red) line 
stands for n = 5 x lO^^cm"^, the dashed (green) line for 
n = 2.5x 10"cm"^ and the full (black) line for n = lO^cm"^ 



given hy m — n but by a lower value. Nevertheless, at 
low densities (n < lO^^cm"^) the picture is not modified 
by self-energy corrections and we expect ferromagnetic 
ordering at T = 0. 



In this Section, we calculate the bilayer density-density 
correlation function, i.e., the response of the system to an 
external potential. The bare polarization as a function 
of momentum q and frequency oj is given in terms of the 
bare particle-hole bubble shown in Fig. 4 (b). It reads: 



P«(q,c.) = 



(2^) 



nF{E{k))-nF{E{k + q)) 



E{k)-E{k + q) 



i5' 



where np{E) is the Fermi-Dirac function. Hence, the 
imaginary part of the retarded response function is given 
by: 



IjnP(i)'^°t(q,t^) = — / SknF{E{k)) 
2vr j 



(28) 



xY,^5{ioj-{E{k)~E{k + q))) 



Assuming a small electronic density and hence elec- 
trons with momentum k such that |fc — fcol ^ ^Oi the 
energy dispersion can be approximated as: 



E{k) 



\k\ 



1) 



(29) 



where we have set fcp = 1 (thus we measure momentum 
in units of fco) and take afcg as our unit of energy (only 
in this section). The polarization is an odd function of 
Lj and hence we can consider only the case of > 0. 
Defining the function: 

f^{x)=^u-{E{k)-E{k + q)) 

= £,uj + q^ + 2k + 2kqx - 2y^k'^ + q'^ + 2kqx , 



with = ±1 and its zeros: 
« ^ 2kq 



V-Cc^ + (fc - 1)2 ± 1 -k^-q 



the imaginary part of the polarization can be written as: 

r-l+h 



ImP(i)^"*(q,u;) = |-^ ^ dfcfc ^ ^ 



(30) 



?,7=± 
1 



where 6 <C 1 defines the Fermi sea, i.e., b = (fcmax — 
^min)/2 and no explicit spin-dependency is considered. 
The integral in (30) can be evaluated analytically but 
the expressions are lengthy. Here we focus on the limits 
of high- and low-frequency relative to the Fermi-energy, 
ep ~ b'^. In the low- frequency limit, we further consider 
the cases of jql ^ 1 and |q — 2| <C 1, representing forward 
and backward scattering processes with small momentum 
transfer, respectively. 
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For the analytical approximations, it is convenient to where the upper and lower bounds in the integral are 
perform separate substitutions for ^ = +1 and ^ = — 1. given by +(po and —ipo , respectively, where (fo = 
For ^ = —1, we set: sinh~^ hj^fu). We thus obtain 



(fc — 1) = sinh(/3 , 



(31) 



Imp(i),ret,€=-i(q^^-) = J difY^ \/^cosh<^(\^^sinh<^ + 1) 



J y/uj cosh Lp\ 



■7V^cosh<^|| yJDxDi ' 



(32) 



where we abbreviate the denominators as, 
D\ = — (\/a;sinh(^ + 1 — g)^ + (1 + 7\/a;cosh(p)^ , 

£>2 = (x^^sinh(^ + l + g)2-(l+7^cosh(p)2 . |fc - 1| = V^cosh^ . (33) 

The final result must be real in such a way that there 

might be a constraint in the integration domain, denoted The upper integration limit is given by <^o = 
by D. 



For ^ = 1, we only have a contribution if a; < (A;— 1)^ < 
6^ and we can perform the substitution, 

|fc — 1| = \/a; coshi^ . 

gration limi 
cosh~^ b/^/uj. We thus obtain: 



Inip(i),ret,c=i( )^ 2 j y ^si^hy)y(CV^cosh<p + l), |l+7V^;^nhc^| 1 

' 27r^.^>o "^^ti ^|l-|l+7\/^sinh(^|| V^^' ^ ' 



r 



where we abbreviate the denominators as, 

Di = — (^-v/wcoshi^ + 1 — g)^ + (1 + 7\/a;sinh<^)^ , 
D2 = (C\/w coshi^ + 1 + g)^ — (1 + 7\/a;sinh<y?)^ . 

A. High frequency 

At high frequencies, ui ^ b^, we can expand sinhc/j ^ 
(f, cosh(p ~ 1 and the resulting integral can be per- 



formed. Both denominators, Di and D2, have to be real, 
altering the integration domain. But for 6^0, the lower 
and upper bounds are given by ^.bj \fuj or zero (since the 
crossover region is missing). In this approximation, the 
integrand can thus be expanded in ip and only the con- 
stant term can be kept. We obtain: 



J 



ImP(^).-(q,.) = -(26) A ^ l^ + ^^L. Re 

7=± ' ' 



1 



=Re- 



1 



+ 7V^II V(l + 9)'-(l+7V^)' V-(l-9)2 + (l+7V^)' ■ 
I 



(35) 



The restrictions, i.e., the regions in which ImP(i)'^«^* = 0, For w ^ 0, we obtain: 

are obtained by imposing real denominators. The results 

are shown on the left hand side of Fig. 8. ImP^ ^^''^ (q,a;) ^ f io-^l'^l{q^\ - {q/2ko)'^) ; < g < 26 

-(7r/26fco) I ; 9 > 26 

(36) 

We thus observe a pronounced non- Fermi liquid behavior 
in the high-energy regime, b^ < w < k^. Only in the small 
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-5 Im P +0. -0.5 Im P +0. 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

q / ko q/b 



FIG. 8: Left hand side: Density plot of the imaginary part 
of the polarization, ImP''^'''°'(q, as function of q and lj 

in units of ka for 6^0. Right hand side: Density plot of 
the imaginary part of the susceptibility, ImP'^^''^'^* (q, a;), as 
function of q and uj in units of 6 ^ fco- 



-47rgImP(i)'"*(q,w)): 

r ; < g < 2& 

P^l (c./62)i/2/V2 ;<z = 26 (39) 

[ {lo/9){1 + l/^l-{2b/qY) ■q>2h 

The results are shown on the right hand side of Fig. 8. 
One sees a linear mode for small wave vectors q 1, 
reminiscent of ID Luttinger liquids. Nevertheless, the 
asymptotic behavior of the polarization for small energies 
shows Fermi-liquid behavior for q ^ 26. At q = 26, we 
find non-Fermi liquid behavior since two parallel regions 
of the Fermi surface are connected. The non-analyticity 
at q = 26 leads to Friedel oscillations with period n/h 
decaying as at large distances. These results are con- 
sistent with alternative analytical treatments for generic 
Fermi surfaces with arbitrary curvature^'^'^^. 



energy regime uj <b^ , Fermi-liquid behavior is recovered 
(see next section). 



B. Low frequency; forward scattering 

In the low- frequency regime, a; ^ 6^, we limit ourselves 
to forward scattering processes with q<^l. Only consid- 
ering the lowest order of ^/ui ^ q, we obtain for <^ = — 1 
the following approximate denominators: 



Di^ 2q + ^2y/uje 



D2->2q- 72^/^e 



yielding, 

Imp(i).rot,«=-i(q^^) 

/'min(ln g^jSinh"" "'" b) 
- sinh~ -'- b 



(37) 



2n 



dip 



2ip 



where we defined b — b/ ^fu) and q — qj ^/lo- 

For ^ = 1, and again only considering the lowest order 
of ^ g, we obtain the following approximate denom- 
inators: 

Di^2q- C2^/ZJe-''^'^ , D2^2q + Q2^e-^^'f , 
leading to, 

Inip(i).r<=t'«=i(q,u;) (38) 

^m.in(ln g,cosh~ "'" &) 



2^ 7o 



1 



1 

2^ 



cosh b 



1 



r,2 



max(— In g.O) 



The analytical expressions are presented in appendix 
B. The behavior for ^ is given by (P = 



C. Low frequency; backward scattering 

For q « 2ko, i.e., for backward scattering processes, 
we consider scattering processes with momentum 2 -I- g 
and \q\ <C 1. We then get for ^ = — 1 the following 
approximate denominators for low energies ^/w ^ q: 



For <^ = 1, we obtain, 



Di T2q + C2y^i 



Do 



Do 



(40) 
(41) 



The analytical solution is presented in appendix B. 
Defining + = -47rImP(i)'™*(2k^, w) with |2kJ| = 
2 + q, where q can be positive or negative, we obtain the 
following behavior for uj 0: 



1 



1 



2^ 



(42) 



(c./62)/2^^^ 

{Lu/b^y/yV2b 

(c^/62)/2/M 

(c^/62)i/2/V26_ 

(c^/62)/2V26^ 

(6^/62)1/2/726 





q < -26 

q = -26 

-26 < g < 

q^O 

< g < 26 

g = 26 

q>2b 



(43) 



The result for "positive" and "negative" back- 
scattering is shown on the right and left hand side of 
Fig. 9, respectively. 

There are three wave numbers that connect two par- 
allel regions of the Fermi surface and thus lead to non- 
Fermi liquid behavior. This results in Friedel oscillations 
with period Tr/feg modulated by oscillations with period 
7r/6, both decaying as at large distances. These are 
enhanced by a factor (26)~i/2 compared to the Friedel 
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-0.5 Im P +0. 



-0.5 Im P +0. 



o/b2, 




-0.5 0.5 1 1.5 2 2.5 3 



q/b 



q/b 



FIG. 9: Density plot of the imaginary part of the polarization, 
Ijjip(i),rct(^2kJ,tj) with |2kJ| = 2 + g, as function of q and uj 
in units of 6 ^ fco. 



Numerical 




Numerical 


— Analytical* 1.5 




— Analyiical/1.3 




FIG. 10: (Color online) The imaginary part of the polariza- 
tion, ImP<^''"'=*(<3,c^), for q = 36/2,36 (left hand side) and 
q = 2ko — b/2, 2ko — 3b (right hand side), as function of uj. 
The solid (black) line shows the numerical and the dashed 
(red) line the analytical result multiplied by a global factor. 



oscillations originating from forward scattering processes 
with q = 2b. 

We have integrated Eq. (30) numerically and obtain 
good agreement between the numerical and the analyti- 
cal result of appendix B in the whole energy regime and 
for all wave numbers up to a global factor 0.75 and 1.3 
in case of forward and backward scattering processes, re- 
spectively. In Fig. 10, the numerical (full line) and ana- 
lytical (dashed line) results for ImP^^)'™* are shown for 
q = 3b/2, 3b (left hand side) and q = 2fco - 6/2, 2ko - 36 
(right hand side) as function of with uj. In all cases, 
the low-energy regime is characterized by Fermi-liquid 
behavior whereas the high-energy shows the square-root 
divergence predicted by Eq. (36). 



D. The plasmon spectrum 

Within the random phase approximation (RPA), the 
electronic dielectric constant is given by: 



6(q,c.) = l-i,qF(i)'-'(q,^), (44) 

2/1^1 T^lic plasmon dispersion, f2g, 
is then given by the zeros of the dielectric function, 
e(q, rig) = 0. In the long wavelength limit, q — ^ 0, we 
have: 



where Vq — 27re'/|q 



2q^an 2q'^Vvpn 



2+2 



(45) 



Let us, at this point, recall the result for a 2D electron 
gas^^ where the dispersion is given by E{k) = fc^/ {2m*), 
where m* is the effective mass. The polarization function 
is given by: 



^'iDEG (q^^) 



nq 



leading to the well-known result: 



/ 2W ^ 



(46) 



(47) 



Coming back to the Fermi ring and introducing the ef- 
fective mass m* — (Aa)^^ , we obtain the same expression 
if we include the degeneracy factor gv = 2 for the two 
nonequivalent Fermi points in one graphene layer. With 
the electrostatic relation V = e^nd {d being the distance 
between the two layers), we thus obtain 



= "^^ '"'^ ^/dq 



V. CONCLUSIONS 



(48) 



We have analyzed the effect of electron-electron inter- 
actions in a biased bilayer graphene in the regime where 
the Fermi surface has the shape of a ring. We have stud- 
ied the stability of this system towards ferromagnetic or- 
der using a one-band model of the "Mexican hat" dis- 
persion. We find that the spin polarized phase is stable. 
Unfortunately, our single band calculation does not allow 
the calculation of the saturation magnetization, which we 
plan to study in the future, taking into account the full 
band structure. 

The unusual electronic occupation in fc-space also leads 
to a deviations from the predictions of Landau's theory of 
a Fermi liquid, ImP ~ |a;| at low energies, {Vvp/t']^)b'^ <C 
e < {Vvj,/tl)k^, (see Fig.[10]). The ImP - |e|-i/2 jg. 
pendence found in this range translates into a quasipar- 
ticle lifetime which decays as r(e) cx ^/\e — efP^- 

At low energies, e <C {Vvp/t\)b'^, the system resem- 
bles a Fermi liquid except for wave numbers which con- 
nect two parallel regions of the Fermi surface, i.e., for 
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q = 2b, 2ko and q = 2fco ± 2b. Note that, as the Fermi 
velocity is proportional to the width of the Fermi ring, 
6, the imaginary part of the response function grows as 
b~^. The existence of two Fermi lines implies Friedel 
oscillations with period w/b. The plasmon dispersion 
shows typical features of two-dimensional systems, nev- 
ertheless the energy scale is greatly enhanced compared 
to a 2DEG. 



APPENDIX B: ANALYTICAL EXPRESSION FOR 
THE SUSCEPTIBILITY 

Wc will here present the approximate analytical result 
of the imaginary part for the polarization. For this we 
define the following abbreviations: q = q/y/tJ, b = bj \fili, 
hs = h+ \/\P- + 1 and 6c = 6 + \/\P- — 1 using sinh"""^ x = 
\n{x + Va;^ + 1) and cosh~^ x = \n{x + \/x^ — 1). 
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APPENDIX A: COULOMB INTEGRALS 

In this appendix, we list the results for the two relevant 
Coulomb integrals. The self-energy corrections involve 
the following Coulomb-integral: 



h{a,y) 



27r <.! 



xdxdO 



(Al) 



J a \/ x"^ + y"^ zL 2xy cos 
= 2{l + y)E[g{y)]+2{l-y)K[g{y)] 
- 2{a + y)E[f{y, a)] - 2{a - y)K[f{y, a)] , 

with g{y) = 4j//(l + yf and /(y, a) = Aay/{a + yf. 

The exchange energy gives rise to the following 
Coulomb- integral: 



4 r 

a) = J ydyh{a,y) = ^^^-^[2{l + a){l + a 



3(1 + a) 

- (1 + af{l + a^)E[h{a)] + (1 - a^fK[h{a)] ' 

(A2) 

with h{a) = 4a/(l -ha)2. 

Both integrals involve the elliptic integrals E{m) and 
K{m) defined as 



E{m) 



r/2 



{l-msin^ey/^de, (A3) 



and 



K{m) 



I.-K/2 

= / {l-msm^6)-^l'^d6. (A4) 
Ja 



For the analytic representation of the polarization 
of forward-scattering processes, we define f{x) = 
cosh^^(q'x) and abbreviate P = — 27rg'ImP(^)'™*(q, w). 
For < g < 6, we then have 



P 





fibs) - f{bc) 



LO > q{2b + q) 

q{2b-q) < u < q{2b + q) 

< w < q{2b - q) 



For b < q < 2b, we have 



P= < 





fibs) 

fibs)- fibe)+f Fib-') 

I fibs) - fibc) 



uj > qi2b + q) 

b"^ <uj < q{2b + q) 
q{2b-q) <LJ <b^ 
0<w < g(26 - q) 



For q > 2b, we have: 

P = fibs)eiqiq + 26) -co)- fib-')Qiq{q - 26) - u) 
-[/(6,)-/(6-i)] 6(6^-0;) 

2. Backward scattering 

For "exact" back-scattering q = 2ko, we have with 
P = — 27ra;^/^ImP(^)'''''*(2ko, w) the following expression: 

P = Vb^ - Vbj' - \/bl — Vbj' 

For "positive" back-scattering, i.e., q = 2kQ with 
|2kQ I — 2 = q > 0, we define fix) = cos~'i^/q/x) and 
abbreviate P = -2nq'/'^ImP^'^''<'\2l4 ,w). 

For q < b, we have; 

P = fibs)-fib;')eiu;~qi2b + q)) 

- [fibc) - fib-')eiio - g(26 - q)] 9(6^ - a;) 
For b < q < 2b, we have: 



fibs)-fib7') 
P={ fibs) 

fibs) - fibc) 

For q > 26, we have: 

fibs) - fibj') 



P 



fibs) 




u) > qi2b + q) 

g(26 -q) <uj < g(26 + q) 

< w < g(26 - q) 



UJ > qiq + 26) 

qiq - 26) < w < qiq + 26) 

0<ui<qiq- 26) 
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For "negative" back-scattering, i.e., q = 2kQ with 
2 — |2kQ I = q > 0, we define fsix) = sinh~^ {\/qx) 
and fc{x) = cosh~^ (y/qx) . The polarization is given by 
pi + p2 ^ -27rgi/2lmP(i)'^<'*(2ko ). The first part is 
independent of the relative value of q and ui and reads 

= fsibs) - fsib;') - [/.(6c) - Mb-')] 6(62 _ ^) 

For < g < 6, we further have: 

r ; w > g(26 + g) 

= < fc{bs) ; «(26 - g) < < g(26 + q) 

[ /c(6.) - /c(6c) ; 0<oj<q{2b-q) 



For a < q < 2b, we have: 



' ; w > q{2b + q) 

fc{bs) ; 62 < a; < g(26 + g) 

/c(6«) - /c(6c) + /c(6-i) ; g(26 -q)<io<b^ 

[ fcibs) - fc{bc) ■,0<oj<q{2b-q) 



For gr > 26, we have: 



= fc{bsMq{q + 26) -uj)- fc{b;')&{q{q - 26) - to) 
- [fF{bc)-fF{b-')]e{b'-uj) 
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